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Abstract 

We consider diagonal and non-diagonal QCD sum rules for the ground state 
heavy baryons to leading order in 1/rriQ and at next-to-leading order in as- In the 
non-diagonal case we evaluate the eight different two-loop diagrams which determine 
the perturbative as-corrections to the Wilson coefficient of the quark condensate in 
the Operator Product Expansion. The QCD corrections to the non-diagonal sum 
rules are moderate compared to the QCD corrections in the diagonal case. We also 
consider constituent type sum rules using constituent type interpolating currents. 
The obtained results are in reasonable agreement with the corresponding results 
obtained in the diagonal case. As central values for the bound state energies we 
find ?h,(Aq) — toq ~ 760 MeV and m(Sq) — toq ~ 940 MeV . The central values for 
the residues are given by F(A Q ) ~ 0.030 GeV 3 and F(Eq) ~ 0.038 GeV 3 . 
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1 Introduction 



The knowledge of the non-perturbative properties of heavy hadrons such as their binding 
energies or their weak transition matrix elements are of fundamental importance for the 
determination of the fundamental parameters of the Standard Model. Among these are 
the quark masses and the values of the Cabibbo-Kobayashi-Maskawa matrix elements. 
A convenient access to the properties of heavy hadrons containing one heavy quark is 
given by the Heavy Quark Effective Theory (HQET) which provides a systematic power 
series expansion of physical matrix elements involving heavy hadrons in terms of the 
inverse of the heavy quark mass (see for example [I]]). While the case of heavy meson 
systems has been analyzed in great detail, corresponding calculations for heavy baryon 
systems have been lagging behind. This is unfortunate since results from the analysis 
of heavy baryon systems are expected to provide important supplementary information 
on the non-perturbative dynamics of QCD and on the fundamental parameters of the 
Standard Model. The importance of further theoretical studies on heavy baryon systems 
is highlighted by the fact that there is now an abundance of new experimental data on 
heavy baryon decays sparked by recent advances in microvertexing techniques. This data 
needs to be analyzed and interpreted theoretically. 

A convenient and well-trusted tool to investigate the non-perturbative properties of 
heavy hadrons is the QCD sum rule method @. The first application of the QCD sum 
rule method to heavy baryons was considered some time ago by Shuryak || who studied 
heavy baryons in the static limit given by the leading term in the 1/rriQ expansion. The 
work of Shuryak || was revised and extended in ||, [|. An analysis of heavy baryons 
containing large but finite quark masses ttiq was undertaken in 11,0, §] . 

We have recently been embarking on a program to improve on previous analysis of 
heavy baryon sum rules by including first order radiative corrections in the analysis. In || 
we determined the two-loop anomalous dimensions of the static heavy baryon currents. 



In [10 1 we determined the perturbative as corrections to the leading, dimension zero, term 
in the operator product expansion (OPE) of the static heavy baryon current correlator. 
Similar to the heavy meson case investigated e.g. in O, the radiative corrections to 
the perturbative dimension zero term are quite large. The results of || [K| were used to 
construct and analyze so-called diagonal QCD sum rules for heavy baryons fllPf . 

Here the term "diagonal" refers to a particular feature of heavy baryon currents and 
their current products. For every baryonic state there are two independent interpolating 
currents even in the static limit |}| [|, ||. One can thus construct diagonal sum rules 
from current correlators of the same baryon current and non-diagonal sum rules from 
current correlators of different baryon currents. The structure of the sum rules for the 
two cases is qualitatively quite different. Nevertheless, they must be considered on the 



same footing. In JT0[ we provided a detailed analysis of the diagonal sum rules. The main 
part of the present paper is devoted to an analysis of the non-diagonal sum rules including 
radiative corrections. We compare our results with those obtained from the analysis of 
the diagonal sum rules. Using results from [[TIJ we also analyze mixed sum rules where 



we use constituent type current combinations in the current correlators. 

In order to provide a brief synapsis of the structural differences of the diagonal and 
non-diagonal sum rules let us briefly recapitulate the main features of the diagonal sum 
rule analysis §, g 0: 

• QCD sum rules based on diagonal correlators feature a leading order spectral density 
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which grows rapidly with energy as p(u) ~ u 5 . This rapid growth introduces a 
strong dependence of the results on the assumed value of the continuum threshold. 



• The QCD radiative correction to the leading order spectral density amounts to 
about 100% at the renormalization scale /i = 1 GeV. 

• Despite of the large QCD corrections to the Wilson coefficients in the OPE, the 
lowest order sum rules and the radiatively corrected sum rules predict nearly the 
same values for the masses and the residues, while the stability region of the sum 
rule results appears at slightly shifted values of the continuum threshold. 

It is clear that one also needs to analyze the non-diagonal sum rules in addition to 
the diagonal sum rules if only for reasons of consistency. A welcome property of the non- 
diagonal correlator is the quite "normal" behaviour p(u) ~ u 2 of the spectral density and 
the fact that the QCD corrections are moderate. 

The paper is organized as follows. In Sect. 2 we introduce our notation and construct 
the correlator of two heavy baryon currents. We also recall the form of the known QCD 
corrections to the dimension zero term of the OPE. In Sect. 3 we present our results on the 
QCD corrections to the dimension three contribution in the OPE, which is proportional 
to the vacuum expectation value of the product of the quark and the antiquark field. We 
also construct generalized QCD sum rules which incorporate both the diagonal and the 
non-diagonal case. Sect. 4 contains the results of our numerical analysis, our final numbers 
and a discussion of the results. In an Appendix we collect our results on the evaluation of 
the radiative two-loop corrections to the dimension-three condensate contribution. The 
results are given for D-dimensional space-time using the most general baryon current 
structure. 

2 Correlator of two baryonic currents 
2.1 Basic notions 

In this section we briefly recapitulate the basic notions involved in the construction of 
QCD sum rules for heavy baryons. This also serves to introduce our notation which 
closely follows the one used in || |K|. The starting point is given by the correlator of two 
baryonic currents = 1,2) 



where k^ is the residual momentum of the heavy quark and is the four-velocity of the 
heavy baryon, the product of both being denoted by uj. The residual momentum and 
the four- velocity are related by p M = rriQV^ + k^, where denotes the momentum of the 
heavy quark and tuq is its mass. As was mentioned before, there are two possible choices 
of interpolating currents for each of the heavy baryon states. Neglecting the flavour and 
colour structure for the moment, these are given by0 



4 Here we use a rather symbolic notation. The Dirac strings T and V can carry Lorentz indices. A 
contraction on the Lorentz indices is always implied when writing the currents in the form of Eq. (0) 




(1) 



J x = [q T Cr iq ]r'Q and J 2 = [q T CT 2 q}T'Q, 
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Table 1: Specific values of the parameter pair (n,s) for particular cases of the light-side 
Dirac structure T. 



where 

Ti = T and T 2 = I>. (3) 

To be more precise, it is clear that the multiplication of the light-side Dirac structure T 
with t does not change the quantum numbers of the interpolating current but does change 
the structure of the current. Multiplying the heavy-side structure V with t, however, does 
not change the structure of the interpolating current since i)Q = Q in the static limit. 

In the static limit one has two types of heavy ground state baryons depending on 
whether the light diquark system is in a spin or in a spin 1 state. We shall employ a 
generic notation and refer to the first type (spin diquark) as Ag-type heavy baryons. 
The Dirac structure of the interpolating current is given by T = 75 and r" = 1 in this case. 
In the second case (spin 1 diquark) one has a doublet of degenerate Sg-type states with 
overall spin 1/2 and 3/2. For the spin 1/2 Eg-type state the interpolating current is given 
by T = 7^ = 7^ — tv^ and V = 7^75. The explicit form of the spin 3/2 interpolating 
current (Sg-type state) can be found e.g. in |L0] . 



For a general analysis it proves convenient to represent the general light-side Dirac 
structure of the currents in Eq. (|^) by an antisymmetrized product of n Dirac matrices 
p = 7^ • • • 7^]. When calculating the one- and two-loop vertex corrections to the baryon 
currents in Eq. (^|) one encounters 7-contractions of the form 7^7" and i)Ti). The j a - 
contraction leads to an n-dependence according to 

la rr = {-l) n {D - 2n)T, (4) 

where D denotes the space-time dimension. The ^-contraction depends in addition on 
the parameter s which takes the value (s = +1) and (s = —1) for an even or odd number 
of #'s in T, respectively. The -^-contraction reads 

ijTi, = (-l) n sr. (5) 

Some of our results in the next sections are given in terms of the most general Dirac 
structure of heavy baryon currents involving the parameters n and s whose definitions 
should be kept in mind. For the convenience of the reader we list the relevant (n, s)-values 
for the cases studied in this paper in Table |I[ 

2.2 Anomalous dimensions of heavy baryon currents 

The one-loop and two-loop renormalization of the static heavy baryon currents and their 
anomalous dimensions were considered in ||] and 0, respectively. Note that the anoma- 
lous dimensions of the currents differ in general from those in conventional QCD. We 
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define the anomalous dimensions in terms of the expansion 7 = Y^k{ a S I ^Tlk- The one- 
loop anomalous dimension for the general current case depends only on the parameter n 
which specifies the light-side Dirac structure. The one-loop anomalous dimension is given 

by |, I 

7i = -^-2) 2 + 2). (6) 

The general (n, s)-dependent formula for the two-loop anomalous dimension case is rather 
lengthy and will therefore not be listed here. The general formula can be found in ||. 
Here we specify to the case of the heavy ground state baryons and give the expansion to 
two- loop order using the MS-scheme and a naively anticommuting 75. One has 
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«15.81 

where the numerical values are given for the case of three light flavours (Nf = 3). 



3 Diagonal, non-diagonal and mixed correlators 

As mentioned before, the two independent currents give rise to two independent types 
of correlators, namely the diagonal correlators {J\J\) and (J2J2), and the non-diagonal 
correlators (J1J2) and (J2J1). In the general case, one may even consider correlators built 
from a linear combination J = aJi + (1 — a) J 2 of currents with an arbitrary coefficient 
a (0 < a < 1). We shall, however, not discuss the most general linear combination 
of currents in this paper. Later on we investigate the case a = 1/2. The choice J = 
(Ji + J2V2 corresponds to a constituent quark model current which has maximal overlap 
with the ground state baryons in the constituent quark model picture. 

Following the standard QCD sum rule method f2|, the correlator is calculated in the 
Euclidean region —uo ?a 1 — 2 GeV including perturbative and non-perturbative contribu- 
tions. In the Euclidean region the non-perturbative contributions are expected to form a 
convergent series. The non-perturbative effects are taken into account by employing an 
OPE for the time-ordered product of the currents in Eq. ([!]). One then has 

(T{J(x),J(0)}) = Y,C d {x 2 )O d 

d 

= C (x 2 )O + C 3 (x 2 )O 3 + C,(x 2 )O i + a o (x 2 )O 5 + ... (11) 

where the are vacuum expectation values of local operators whose dimensions are 
labelled by their subscripts d. O = 1 corresponds to the so called perturbative term, 
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0% = (qq) is a quark condensate term, O4 = ots(G 2 ) is a gluon condensate term, O5 = 
gsiq^nuG^q) is a mixed quark-gluon condensate, and so on. The expansion coefficients 
Cd{x 2 ) are the associated coefficient functions or Wilson coefficients of the OPE. 

A straightforward dimensional analysis shows that the OPE of the diagonal correla- 
tors (J\J\) and (J2J2) contains only even-dimensional terms, while the OPE of the non- 
diagonal correlators contains only odd- dimensional terms. This classification is preserved 
when radiative corrections are included, assuming the light quarks to be massless. We 
apply radiative corrections only to the leading terms in the OPE because the non-leading 
contributions are small. 

The diagonal case was studied in detail in ||10|| . It was shown that the QCD corrections 
to the spectral density of the correlator function P(u) are quite large. It is quite intriguing 
that the contributions of the four different three-loop diagrams that contribute to the 
perturbative dimension zero piece shown in Fig. 1 can be collected into one compact 
formula |10| : 



p^ CD (^,/i) = as 



In ( £-) -(n 2 - 4n + 6) + — (60C(2) + 38n 2 - 137n + 273) 



.(12) 



The number n specifies the light-side Dirac structure of the baryon currents as before. 
Note that the coefficient of the logarithmic term coincides with the one-loop anomalous 
dimension of the diagonal correlator which in this case is equal to two times the anomalous 
dimension of the baryon current itself. 



3.1 Non-diagonal correlators 

The non-diagonal correlator of the two heavy baryon currents reads 

n 12 (u;) = i J d*xj kx (u\T{j x {x)j 2 {vi)}\V) = r;i±lr' 2 ~ tt{y x y 2 )p X2 {u). (13) 

We have suppressed the flavour and colour labels in Eq. ([13]). The OPE for the non- 
diagonal correlator contains a term O^fi) = (qq) proportional to the quark conden- 
sate, a mixed quark-gluon condensate term 5 (/i) = gsiq&fivG^q) = ml(qq), a term 
7 (/x) = (qq)(asG 2 ), and a term Og = as(qq) 3 - Taking into account these four conden- 
sate contributions, the Fourier transform of the scalar correlator function P(w) is given 

by 

P{t) = Pope(*) (14) 

= ~'^W [° M + - 2)°^ + 288 i 1 - 2> M - 9T2° M ) • 

where c is a Clebsch-Gordan type factor which takes the values c = 1 for the AQ-type 
and c = —1/3 for the S^-type doublet {Sq,Sq} ground state baryons. The correlator 
function Pope(^) satisfies the dispersion relation 

poo piuj'^dui 1 

P ope (lu) = P(uj) = / HK ' + subtraction, (15) 

Jo uj' — uj — 1O 

where p(u) = Im(P(u))/'K is the spectral density of the scalar correlation function. Taking 
into account the above four condensate contributions, the lowest order spectral density of 
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the first two contributions is given by 



(qq) 



-IQ 



7T" 



and Pb(uj) = 2(1 



16tt 2 



(16) 



Next we compute the radiative corrections to the quark condensate term p^(uj). There 
are altogether eight different contributing diagrams which are shown in Fig. 2. Their 
contributions were evaluated with the help of the algorithm developed in [13]. As a check 
on the calculation we used a general covariant gauge for the gluon. The gauge dependence 
was found to drop out in the sum of the contributions. 

Collecting together the one- and two-loop contributions to the dimension three scalar 
correlation function one has 



(qq) 

2tt 2 



-2lu 



D-A 



C D + 



9s 



-2w 



2D-8 8 



(4tt)^/ 2 



P 



(17) 



i=l 



where D = 4 — 2e is the space-time dimension. There are a number of colour factors in 
Eq. (0) the values of which are given by C = N c \, d = C 2 = C 3 = C 4 = C 5 = -N C \C B 
and Cg = C-j = C$ = N c \Cf, where N c is the number of colours, Cf = (iV 2 — l)/2N c , and 
Cb = (N c + 1)/2N C . Explicit forms of the scalar coefficients Di defined in Eq. ([17]) are 
listed in Appendix A. When Eq. fllTD is expanded in terms of a power series in 1/e, one 
obtains 
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2tt 2 
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where we have now substituted explicit values for the colour factors with N c = 3. The 
spectral density p^{oj) is given by the absorptive part of Ps(uj). In renormalizing the 
spectral density ^3(0;) one has to take into account both the renormalization factor of the 
baryon currents and the renormalization factor of the quark condensate, 



ps( U ) = Zj.Zj^pT^) 



(19) 



with 10, 12 



"qq 
Zj 2 



1 + 3 



Aire ' 
1 + ^(n 2 

l + ^((n-2 + S ) 2 + 2). 



An + 6) and 



(20) 
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After multiplication with the Z-factors the leading 1/e-contribution in ps{uj) is cancelled. 
The renormalized spectral density is given by 



where 



with 



P3 Born H 
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1 + S 7 "^/^ 



7T- 



u and r (uj/p) = r\ In ( — I + r 2 



P 



2u 



(21) 



:22) 



4 

-(2n 2 - 8n + 7 + 2(n - 2)s) 
3 



and 



(8n 2 - 28n + 37 + 8ns - 14s + 8((2)). 



(23) 



Note that the coefficient 7*1 of the logarithmic term in Eq. ( p2]) coincides with the sum of 
the one-loop anomalous dimension of J\ and J 2 minus the anomalous dimension of the 
quark condensate. The reason is that the same coefficient is involved in the cancellations 
of the 1/e-pole in Eq. flU). 

Explicit values for the correction to the spectral density for the cases of the Aq- and 
SQ-type ground state baryons are given by 



r A (u/p) = 4 In 



and 



r s (w//i) 



In 



A. 

2lv 

jj_ 

2u 



+ -(23 + 8C(2)) 



;(n+8c(2)). 



(24) 



(25) 



The radiative corrections can be seen to amount to about 40 — 60% at the renormalization 
scale \i — 1 GeV. Due to the hermiticity of the current correlator ILj the coefficients r\ 
and T2 do not depend on which of the two non-diagonal current products J\ J2 or J 2 ^i ar e 
taken. 



3.2 Non-diagonal sum rules 

As usual we construct QCD sum rules by invoking parton-hadron duality, i.e. we equate 
the theoretical contribution to the scalar correlation function P(u) given in Eq. (|13D with 
the dispersion integral over the contributions of hadron states. These consist of the lowest 
lying ground state with bound state energy A plus the excited states and the continuum 
contributions. To leading order in I/trq the bound state energy of the ground state is 
defined by 

^baryon = m Q + A, (26) 

where rriQ is the pole mass of the heavy quark. We assume that the contribution of 
the excited states and the continuum contribution sets in above some effective threshold 
energy Eq and can be approximated by the OPE expression 0. For the hadron-side 
contribution p H s to the spectral density we thus write 

PHsM = PGsM + PcontM, (27) 
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where the contribution of the lowest-lying ground state baryon is denoted by pes and is 
given by 

Pgs (uj) = ^F 1 F 2 5(uj-A). (28) 



The residues Fj (i — 1, 2) appearing in Eq. ( |28| ) are defined by the matrix elements of the 
heavy baryon currents according to 



(0\Ji\A Q ) = F iA u, (0| Jj|Sg) = FiY,u and (0\J?\Z* Q ) = -^F^u", (29) 



where u and u v are the usual spin 1/2 and spin 3/2 spinors. Note that F^* coincides 
with FiY, to lowest order of the heavy quark mass expansion which we are working in. 

As is usual we assume hadron-parton duality for the contributions of excited states 
and the continuum contributions. As mentioned before we subsum these contributions by 
defining an effective energy threshold Ec and write p C ont(^) = 0(ui — Ec)p(oj), where p is 
the result of the OPE calculations given in Eqs. (|i~6|) and (|2l|) . With these assumptions 
we arrive at the sum rule 

Ww)= jM-.r^ (so, 

A — uj — l\) JE C uj' — UJ — 20 

or 

A — c<j — 2O Jo — w — «0 
The polynomial contribution Ppc(uj) is defined as the Fourier transform of that part of 
the correlator function P(t) which contains non-negative powers (t 2 ) n (n > 0). Finally 
we apply the Borel transformation 




Bt = lim— — — - ( — — ] n, —uj ^00 (T = —uj/n fixed) (32) 



to the sum rule in Eq. (|3lD . Using Bt(1/ (uj — uj')) = exp(—uj'/T)/T we obtain the Borel 
sum rule 

1 - r E c 

-FMF 2 (p)e- A / T = / p(uj > ,p)e-"'/ T duj' + BP PC (T) =: K(E c ,T } p), (33) 

2 Jo 

where we have reintroduced the /i-dependence of the spectral density which in turn gives 
rise to a p-dependence of the residue. The Borel-transformed polynomial contribution 
BP PC (T) can be obtained directly from Ppc(^) by the substitution t — > —i/T (see the 
discussion in ||). Note that the bound state energy A can be obtained from the sum 
rule in Eq. (|33|) by taking the logarithmic derivative with respect to the inverse Borel 
parameter according to 

dln(K(E c ,T,p)) 

dT- 1 " 1 ' 

Before turning to the numerical analysis of the non-diagonal sum rules we want to 
briefly comment on the scale dependence of the residues. The numerical values given below 
are taken at the specific normalization point p = 1 GeV, while the general dependence 
on the scale p is controlled by the renormalization group equation. At one-loop order the 
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product of the residues at one renormalization point jj 2 can be expressed by the product 
at another renormalization point fi± via 



with 



^1(^2)^2(^2) = F 1 (n 1 )F 2 (ni)U(n 1 ,n 2 ) 



71 /ft 



U(fii,fi 2 ) 



(35) 



(36) 



where U([Ax, fi 2 ) is the evolution function which takes one from the scale fi 2 to the scale 
The coefficients (3\ = 11 — 2/3N c and 71 are the usual first-order terms in the expansion of 
the QCD /^-function and the anomalous dimension ^ P = 7j x +7j 2 — 7?g °f ^ ne non-diagonal 
scalar correlator function P. The two-loop extension of Eq. (p6[) is given by 



exp 



/•as 0*i ) da 7(a) 
lasto) a (3(a) 



71 /ft 



1 + 



4vr A 




(37) 



where /3 n is the n-th order term in the /3-function expansion and 7 n denotes the anomalous 
dimension of the non-diagonal scalar correlator function P at n-th loop order. The two- 
loop evolution function in Eq. Q3"7| ) is obtained as a solution to the renormalization group 
equation including next-to-leading order perturbative terms in as (see also the discussion 
in0§p,0). 

It is evident that we can only extract the value of the product of residues FxF 2 from 
our sum rule analysis. In order to make further progress, we adopt the working hypothesis 
that the residues of the two current options in each case are equal. This assumption is 
corraborated by the results of the diagonal sum rule analysis . This means we replace 
F\F 2 by F 2 in the above formulae when performing the numerical analysis. We note, 
however, that the currents J% and J 2 have different anomalous dimensions and therefore 
F\ and F 2 do not coincide at another renormalization scale fi 2 even if they coincide at the 
scale fix- Returning to the sum rule in Eq. ([33]), one then has 
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Air 



In (Sf)M x c) ~ ffcc) ) rx + f2(x c )r 2 
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where the (n, s)-dependent coefficient functions r\ and r 2 are defined in Eq. 
functions /„ and f l n are given by 



(38) 
|) and the 
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(39) 



In order to simplify the notation we have introduced the abbreviations 



x c ■-- 



Ec F 



^. (E Q f 



7T" 



2N, 



(qq) and (E G ) 



na s (G 2 ) 
32N C (N C - 1) • 



(40) 
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3.3 Mixed sum rules of constituent type 



As mentioned before, we shall not investigate the most general case of mixed sum rules 
but specify to the linear combination of currents J = (J% + J2V2 in the sum rules. The 
light-side Dirac structure of the currents can then seen to appear in the form ~(1 + $)r, 
i.e. one has the projector factor P + = (1 + f)/2 which projects on the large components 
of the light quark fields. In the rest system of the heavy baryon, where = (1; 0, 0, 0), 
this is manifest since then P + = (go) - ^ e Te ^ eT t° this particular linear combination of 
currents as the constituent type current. This linear combination of currents is expected 
to have maximum overlap with the heavy ground state baryons in the constituent quark 
model, i.e. where the light diquark state in the heavy baryon is taken to be composed of 
on-shell light quarks. We mention that the constituent quark model picture emerges in 
the large iV c -limit ||16|| . The tools needed for the sum rule analysis of constituent type 



heavy baryons have been assembled in this paper and in || . The results of the constituent 
analysis are presented in the next section together with the results of the analysis of the 
diagonal and non-diagonal sum rules. 



4 Numerical analysis 

Having the neccessary formulae at hand we next describe our numerical analysis of the 
sum rules and specify our choice of the relevant input parameters. We use the following 



numerical input values for the condensate contributions B [XT 



(qq) = -(0.23 GeV) 3 (quark condensate), (41) 
as(G 2 ) = 0.04 GeV 4 (gluon condensate), and 
gs{q< 7 nvG fll 'q) = rnl(qq) with wig = 0.8 GeV 2 (mixed quark-gluon condensate). 

There are in general two strategies for the numerical analysis of the QCD sum rules. The 
first strategy fixes the bound state energy A from the outset by choosing a specific value 
for the pole mass of the heavy quark and then extracts a value for the residue F. In order 
to obtain information from the sum rules which is independent of specific input values, 
we adopt a second strategy, namely to determine both A and F by finding simultaneous 
stability values for them with respect to the Borel parameter T. 

The first step in carrying out the numerical analysis of the sum rules is to find a sum 
rule "window" for the allowed values of the Borel parameter T. The parameter range of 
T is constrained by two different physical requirements. The first is that the convergence 
of the OPE expansion must be secured. We therefore demand that the subleading term 
in the OPE does not contribute more than 30% of the leading order term. This gives a 
lower limit for the Borel parameter. The upper limit is determined by the requirement 
that the contributions from the excited states plus the physical continuum (even after 
Borel transformation) should not exceed the bound state contribution. This requirement 
is neccessary in order to guarantee that the sum rules are as independent as possible of the 
model-dependent assumptions concerning the profile of the theoretical spectral density, 
i.e. the model of the continuum. 

The lower limit of Ec is given by the requirement that the indicated window should 
be kept open. For the rest, Ec is a free floating variable which is only limited by the 
stability requirements on A and F. 
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Baryon type state 


Ec 


A 


F 


A Q (L.O.) 
Aq (N.L.O.) 

Eg (L.O.) 
Eg (N.L.O.) 


1.2 ±0.1 GeV 
1.1 ±0.1 GeV 
1.4 ±0.1 GeV 

1.3 ±0.1 GeV 


0.77 ±0.05 GeV 
0.77 ±0.05 GeV 
0.96 ±0.05 GeV 
0.94 ±0.05 GeV 


0.022 ± 0.001 GeV 6 
0.027 ±0.002 GeV 3 
0.031 ± 0.002 GeV 3 
0.038 ± 0.003 GeV 3 



Table 2: Results of the diagonal sum rule analysis for the continuum threshold parameter 
E c , the bound state energy A, and the residuum F for Ag-type and Eg-type currents, 
analyzed to leading order (L.O.) as well as next-to-leading order (N.L.O.) 



4.1 Diagonal sum rules 

Let us briefly recapitulate the results of the numerical analysis of the diagonal sum rules 
presented in [|Kj]. The above two requirements limit the allowed range for the Borel 
parameter to 250 MeV < T < 400 MeV. The analysis proceeded in two steps. First we 
analyzed the uncorrected sum rules varying both the continuum threshold and the bound 
state energy. The criterion for the best choice of these two energies is the stability of the 
sum rules with regard to the variation of the Borel parameter T. In the second step we 
included the radiative corrections and again varied both the continuum threshold and the 
bound state energy to obtain the best sum rules stability. The ratio of the continuum 
contribution and ground state contribution depends strongly on the Borel parameter T 
and the continuum threshold energy Ec- Looking e.g. at the sum rule analysis for the 
Ag-type baryons with QCD corrections, the continuum contribution is about 80% of the 
ground state contribution for Ec = 1.1 GeV and T = 250 MeV and then increases with T. 

Because of the new specification for the sum rule window we have repeated the diagonal 
sum rule analysis of [[KJ allowing for slightly different values of A and F. The outcome 
of the numerical analysis is practically unaltered. The values for the Ag-type state can 
be read off from Fig. 3. Fig 3(a) shows the dependence of the bound state energy A on 
the Borel parameter T and Fig. 3(b) shows the dependence of the residue on T, both for 
the leading order sum rule. Fig. 3(c) and Fig. 3(d) show the same dependencies for the 
radiatively corrected sum rules. The same analysis is repeated for the Eg-type states in 
Fig. 4. The results of the numerical analysis both without and with radiative corrections 
are given in Table [| 

4.2 Non-diagonal sum rules 

In the case of the non-diagonal sum rules, the "window" of permissible values for the Borel 
parameter is wider than in the diagonal case and it is given by 250 MeV < T < 600 MeV. 
Proceeding in the same manner as in the case of the diagonal sum rule analysis, we obtain 
best stability values when varying T. The values for the Ag-type state can be read off from 
Fig. 5, and the values for the Eg-type state can be obtained from Fig. 6. For the Ag-type 
baryons the stability appears at values of T and Ec where the continuum contribution 
is about 100%. If we try to decrease this contribution relatively by increasing Ec, the 
stability in T becomes worse. As can be seen from Fig. 5(a), for e.g. Ec = 1.3 GeV 
the contribution of the continuum is less than 40% on the left hand side, but stability is 
lost. These considerations show that the relative error of our estimate can be taken to be 
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Baryon type state 


Ec 


A 


F 


Ag (L.O.) 
Aq (N.L.O.) 

Eg (L.O.) 
Eg (N.L.O.) 


1.0 ± 0.10 GeV 
1.0 ±0.10 GeV 
1.5 ± 0.10 GeV 
1.2 ±0.10 GeV 


0.75 ± 0.10 GeV 
0.72 ± 0.10 GeV 
1.16 ± 0.10 GeV 
0.94 ± 0.10 GeV 


0.024 ± 0.002 GeV 6 
0.032 ± 0.003 GeV 3 
0.045 ± 0.003 GeV 3 
0.039 ± 0.004 GeV 3 



Table 3: Results of the non-diagonal sum rule analysis for the continuum threshold pa- 
rameter Ec, the bound state energy A, and the residuum F for Ag-type and Eg-type 
currents, analyzed to leading order (L.O.) as well as next-to-leading order (N.L.O.) 



Baryon type state 


Ec 


A 


F 


Aq (L.O.) 
A Q (N.L.O.) 

Eq (L.O.) 
E Q (N.L.O.) 


1.1 ± 0.10 GeV 

1.1 ± 0.10 GeV 
1.3 ± 0.10 GeV 

1.2 ± 0.10 GeV 


0.77 ± 0.10 GeV 
0.77 ± 0.10 GeV 
1.03 ± 0.10 GeV 
0.94 ± 0.10 GeV 


0.034 ± 0.004 GeV 6 
0.032 ± 0.004 GeV 3 
0.045 ± 0.004 GeV 3 
0.036 ± 0.004 GeV 3 



Table 4: Results of the constituent type mixed sum rule analysis for the continuum 
threshold parameter E c , the bound state energy A, and the residuum F for Ag-type 
and Eg-type currents, analyzed to leading order (L.O.) as well as next-to-leading order 
(N.L.O.) 



approximately 10%. The situation for the Eg-type baryons is much better. For example, 
for the radiatively corrected sum rules the radio of the continuum and the ground state 
contribution is 50% for the central value E c = 1.2 GeV and 30% for E c = 1.5 GeV at 
the left end of the allowed range for the Borel parameter T. 

The numerical results are given in Table |3|. Assuming relative errors of 10% for the 
bound state energy and 20% for the residue, the obtained values are in agreement with 
the results of the analysis of the diagonal sum rules, where the values for the Eg-type 
baryon are the more reliable one. 

4.3 Constituent type mixed sum rules 

The use of a constituent type interpolating current J — ( J\ + J%)j1 combines the two 
sum rule formulas for the diagonal and the non-diagonal case, taking one half of each 
part. The "window" of permissible values for the Borel parameter T is now given by 
300 MeV < T < 700 MeV. In Fig. 7 we show the results of the sum rule analysis for the 
Ag-type baryons, and in Fig. 8 we show the results for the Eg-type baryons. 

The numerical results of the analysis are given in Table f|. The constituent type sum 
rules show an improved stability on the Borel parameter T as compared to the non- 
diagonal sum rules, but the stability is not as good as in the diagonal case. Within 
the assumed errors the results are again in agreement with both the diagonal and the 
non-diagonal sum rule analysis. 
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4.4 Comparison with experimental values 



Finally we want to compare our results for the bound state energy with the existing 
experimental values for the baryon masses. For such a comparision we need to know 
the pole mass of the heavy quarks which can be extracted from the heavy quarkonium 
and heavy-light mesons [TJ, [L9|, |20[. The quoted value of the bottom quark pole mass 
varies from m b = 4.55 ± 0.05 GeV pf and m b = 4.67 ± 0.10 GeV fH to m b = 4.80 ± 
0.03 GeV 0. Assuming the range 4.6 GeV < m b < 4.9 GeV, the mass m(A 6 ) = 5642 ± 
50MeV^ of the baryon A b [§l] results in a range 740 MeV < A(A 6 ) < 1040 MeV for the 
bound state energy. Our central value A(Aq) = 760 MeV for the bound state energy 
suggests a pole mass of m b = 4880MeV for the bottom quark. 

Taking the experimental results for charm-quark baryons, namely m(A c ) = 2284.9 ± 
0.6 MeV and m(E+) = 2453.5 ± 0.9 MeV our central values A(A Q ) = 760 MeV and 
A(Eq) = 940 MeV predict a mean pole mass of m c = 1520 MeV for the charm quark. 



5 Conclusions 

We have considered the operator product expansion of the correlator of two static heavy 
baryon currents at small Euclidian distances and determined the as radiative corrections 
to the first and second Wilson coefficient in the expansion. Based on the operator product 
expansion we have formulated and analyzed heavy baryon sum rules for the Ag-type and 
Eg-type heavy baryons using two different types of interpolating fields for the baryons in 
each case. In this paper we have constructed and analyzed the non-diagonal sum rules 
built from the correlators of two different currents including radiative corrections. The 
non-diagonal sum rules bring in some new features such as a more "normal" behaviour of 
the spectral density p{u) ~ (qq)^ 2 and moderate QCD corrections to the spectral density 
as compared to the diagonal case. We have taken a second look at the diagonal sum rules. 

We have also set up and analyzed constituent type heavy baryon sum rules where we 
have used interpolating currents that are expected to have a maximum overlap with the 
heavy baryon's light diquark system in the constituent quark model picture. All the three 
types of sum rules show acceptable stability in their dependence on the Borel parameter, 
where the best stability was obtained for the diagonal sum rules. The results of the 
three types of sum rules (diagonal, non-diagonal, constituent type) on the bound state 
energy an the residues of the heavy ground state baryons were found to be consistent with 
each other, where the values obtained for the Eg-type baryons are more reliable than the 
results for the Ag-type baryons. 
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Appendix 



In this appendix we collect our results on the evaluation of the one-loop and two-loop con- 
tributions to the non-diagonal correlators of two heavy baryon currents. The contributing 
diagrams are shown in Fig. 2. Introducing the abbreviation E n = T(l — e) n r(l +ne) (with 
integer numbers n — 1, 2, 3, . . . ) we obtain 

4f £i „ „ 2(D-2)f E 2 /a , 

n _ u 1 r) — n v / u z (A 1 ^ 

" (D-4)(£>-3)' 1 " 2 (L>-4) 2 (L>-3)(2L>-7)' 1 ' 

8(D - 2)f E 2 4(D-2)(3D-10)f ^ 
3 ~" (D-4) 3 (D-3) 2 ( J D-4) 3 ( J D-3) 2 (2 J D-7)' 1 ' 

(D-4)fi + f 2 2(D - 2)f - £>fi + f 2 

Da ~ ( jD -4) 2 ( J D-3)(2 J D-7) 2 ' 5 ~ ( jD -4) 2 ( J D-3)(2D-7) 2 ' (A3J 

2D(D-2)f E 2 2(D - 2)f E 2 

6 ( J D-4) 2 ( J D-3)(2D-7)' 7 ( J D-4) 2 ( J D-3)(2 J D-7)' 1 ' 

D - A{D ~ 2)foE2 (AS\ 

8 ~ (D — 4) 2 (D — 3) 2 (2D — 7) 1 ; 

where r = Tr(f-^r^), Ti = Tr(f7 M r7 M ), and T 2 = Tr(f'07 At 7^r7 iy 7 M ^). 
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Figure Captions 

Fig. 1: Radiative corrections to the diagonal correlator. (0) lowest order two-loop 
contribution, (l)-(4) 0(ots) three- loop contributions. 

Fig. 2: Radiative corrections to the non-diagonal correlator given by the dimension 
three condensate contribution. (0) lowest order one-loop contribution, (1)- 
(8) 0(ot s ) two-loop contributions. 

Fig. 3: Bound state energy and residue of the Aq as functions of the Borel pa- 
rameter T (diagonal case). Plotted are five curves for five different values 
of the threshold energy Eq spaced by 100 MeV around the central value 
Eq = E c cst . Ec increases from bottom to top. 

(a) lowest order sum rule results for the bound state energy A (A) 

(b) lowest order sum rule results for the residue F A 

(c) O(ots) sum rule results for the bound state energy A(A) 
for the current Jm 

(d) 0(as) sum rule results for the residue Fa for the current Jm 

Fig. 4: Bound state energy and residue of the Eq as functions of the Borel pa- 
rameter T (diagonal case). Plotted are five curves for five different values 
of the threshold energy Ec spaced by 100 MeV around the central value 
E c = E c cst . E c increases from bottom to top. 
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(a) lowest order sum rule results for the bound state energy A(E) 

(b) lowest order sum rule results for the residue F s 

(c) 0(0:5) sum rule results for the bound state energy A(E) 
for the current J^i 

(d) 0(as) sum rule results for the residue F s for the current J^i 

Fig. 5: Bound state energy and residue of the Aq as functions of the Borel param- 
eter T (non-diagonal case). Plotted are five curves for five different values 
of the threshold energy E c spaced by 100 MeV around the central value 
E c = E c cst . E c increases from bottom to top. 

(a) lowest order sum rule results for the bound state energy A (A) 

(b) lowest order sum rule results for the residue Fa 

(c) 0(0-5) sum rule results for the bound state energy A(A) 
for the two currents Jai and J A2 

(d) 0(05) sum rule results for the residue F A 
for the two currents Jm and J A2 

Fig. 6: Bound state energy and residue of the Eq as functions of the Borel param- 
eter T (non-diagonal case). Plotted are five curves for five different values 
of the threshold energy Ec spaced by 100 MeV around the central value 
Ec = E c cst - Ec increases from bottom to top. 

(a) lowest order sum rule results for the bound state energy A(E) 

(b) lowest order sum rule results for the residue F s 

(c) 0(05) sum rule results for the bound state energy A(E) 
for the two currents Jsi and J^ 2 

(d) 0(05) sum rule results for the residue F s 
for the two currents Jsi and Js 2 

Fig. 7: Bound state energy and residue of the Aq as functions of the Borel parame- 
ter T (constituent type mixed case). Plotted are five curves for five different 
values of the threshold energy Ec spaced by 100 MeV around the central 
value Ec = E c cst . Ec increases from bottom to top. 

(a) lowest order sum rule results for the bound state energy A(A) 

(b) lowest order sum rule results for the residue Fa 

(c) 0(05) sum rule results for the bound state energy A(A) 
for the two currents Jm and J42 

(d) 0(05) sum rule results for the residue F\ 
for the two currents Jm and Ja2 

Fig. 8: Bound state energy and residue of the Eg as functions of the Borel parame- 
ter T (constituent type mixed case). Plotted are five curves for five different 
values of the threshold energy E c spaced by 100 MeV around the central 
value E c = E c cst . E c increases from bottom to top. 

(a) lowest order sum rule results for the bound state energy A(E) 

(b) lowest order sum rule results for the residue F^ 

(c) 0(05) sum rule results for the bound state energy A(E) 
for the two currents Jsi and J^2 

(d) 0(05) sum rule results for the residue F s 
for the two currents Jsi and Js 2 
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